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We demonstrate that repulsive power law interactions can induce directed transport of particles 
in dissipative, ac-driven periodic potentials, in regimes where the underlying noninteracting system 
exhibits localized oscillations. Contrasting the well-established single particle ratchet mechanism, 
this interaction induced transport is based on the collective behaviour of the interacting particles 
yielding a spatiotemporal nonequilibrium pattern comprising persistent travelling excitations. 
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Introduction The ratchet effect describes the emer¬ 
gence of a directed particle transport due to an interplay 
of symmetry breaking [ 3 ^] and nonlinearity in a system 
where all applied forces an d g radients vanish after aver¬ 


aging over time and space [2|-[7[]. Originally inspired by 
the desire of understanding how molecular motors con¬ 
vert chemical energy into directional motion Q, ratchets 
have now advanced to a widespread and topical paradigm 
with applications in various branches of nonequilibrium 
atomic Riffl , condensed matter BE® and bio¬ 
physics 18 ff Specifically modern atomic cooling and 
trapping techniques 2 J, [25] allow for particularly clean 
and versatile realizations of ratchets in optical lattices 
[jj-fTTl 2(1 27j. Since the established ratchet effect is of 
single particle character that can survive interactions, 
most works focused on the non-interacting single parti¬ 
cle regime (for reviews see liSEl) or on continuous 
models with infinitely many degrees of freedom and ef¬ 
fectively linear or Kuramoto-type interactions [1, 28|, 2(|. 
However, many relevant setups, such as driven dipolar 
atoms or molecules 30j, 31 1 , charged colloids [ 32 I] or mi¬ 
croscopic ions (33^| experience power-law (Coulomb and 
dipole-dipole) interactions, evoking the question of their 
impact on directed currents. For the dilute Hamiltonian 
regime, it was already shown that these interactions can 
induce intriguing dynamical reversals of the transport di¬ 
rection 341. 

In this work, we demonstrate that repulsive power law 
interactions can even be responsible for the emergence of 
a directed particle transport. Considering ions or dipoles 
in dissipative and laterally oscillating lattice potentials, 
interaction induced currents appear in parameter regimes 
where frictional energy losses suppress any directed cur¬ 
rent in the corresponding noninteracting system. The 
physical origin of these currents is an interaction induced 
collective behaviour of the particles which is represented 
by a spatiotemporal nonequilibrium pattern comprising 
persistent travelling excitations and has to be carefully 
distinguished from the well established single particle 


ratchet mechanism. 

Setup We consider N repulsively power law interact¬ 
ing point particles in one dimension with coordinates 
Xj and mass m that are exposed to a frictional force 
Fji = —'yx and a spatially periodic and laterally os¬ 
cillating lattice potential V(x,t) = Vbcos 2 [A:a: — /(f)]. 
Here L = ir/k is the distance between two adjacent min¬ 
ima and /(f) = a [cos(wof) + sin( 2 wof)] is the biharmonic 
driving law, with frequency wo and amplitude a. The dy¬ 
namics of this model is described by N coupled nonlinear, 
non autonomous Newtonian equations of motion: 


N 


mxj + 'yij + Vo/csin {2k [ Xj — /(f)]} — ^ Fjk = 0 

k=l,k^j 

Xj Xk 


Fjk = a t 


| Xj - Xk 


\r+2 


(1) 


This model might describes charged (r = 1) or dipolar 
particles (r = 2 ) in an ac-driven (optical) lattice poten¬ 
tial including a dissipative frictional force. 

To study transport phenomena at a fixed particle density 
we employ periodic boundary conditions x+LN = x. We 
incorporate the corresponding Ewald sums for the inter¬ 
action forces Fjk in our simulations [35], but note that all 
results presented below remain qualitatively valid when 
truncating the interaction beyond next neighbour ones. 
Physically, we can distinguish four different forces de¬ 
termining the dynamics of a given particle ensemble: 
(i) the force exerted by the instantaneous static lattice 
( F s := — Vbfcsin(2/ca;) ~ Vok). If dominant, this force 
drags the particles along the minima of the oscillating 
lattice, {it) The time dependent pseudo force which acts 
on a particle in the comoving coordinate system of the 
lattice (F 0 := -m/(f) -y/(f) ~ (-aw 0 )y/ 7 2 + (mw 0 ) 2 ). 
(in) The frictional forces Ff = —'yx ~ 7 W 0 T slowing 
down the particles in favour to keep them at fixed po¬ 
sitions in the laboratory coordinate system, (iv) The 
repulsive interaction forces Fjk ~ a/L r+1 . In the fol¬ 
lowing parameters are chosen such that F s , F a , Ff and in 
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many cases also Fjk are comparable for sufficiently long 
times, i.e. none of these forces can be treated as a small 
perturbation. To be in a regime where the driving forces 
and the nonlinear components of F s as well as the power 
law shape of Fjk play a significant role, ojq is chosen close 
to the resonance of the underlying linearised single par¬ 
ticle problem (a damped driven harmonic oscillator), i.e. 
2 o>o ~ \J TVF^ j m. This way, we are far from the com¬ 
paratively simple regimes of adiabatic.ally slow and fast 
driving resulting in effectively time independent prob¬ 
lems. 

In the following we investigate the impact of different 
interaction strengths a > 0 on the transport proper¬ 
ties of spatially uniform initial states with Xj(t = 0) = 
L(j — N/2 + 1/2) + a (N even and j = 0, 1...N — 1) and 
random velocities Vj(t = 0), uniformly chosen in an in¬ 
terval (—Vo, Vo). To avoid capture by specific low-energy 
asymptotic states (attractors) we choose vq to be large, 
such that frictional forces dominate at short time scales. 
Experimentally, different values of a can of course be 
achieved indirectly by tuning other parameters such as 
wo, Vo or 7 that are easier to control. 

Results We integrate Eq. (CD using the Runge-Kutta- 
Dormand-Prince integration scheme 36] and calculate 


the time evolution of the ensemble averaged position 
(d(t)) := (1/N) X) Xj(t) (without setting x + NL back 
to x). Although noise is of course naturally present in 
dissipative many particle systems due to the fluctuation- 
dissipation theorem 37, [3||, we do not include it explic¬ 
itly in our simulations but simply note that our results 
are indeed robust with respect to weak (Gaussian) noise. 
Fig. [Hand its inset (a) show the time evolution of ((d(t))), 
where ((.)) refers to the averaging over many JV-particle 
ensembles. For vanishing interactions (a = 0) as well as 
for weak interactions (a = 1.0) (( d(t ))) saturates after a 
short initial growth at t/T ~ 20 and we observe that all 
particles asymptotically perform pure on site oscillations. 
Increasing the strength of the repulsive interactions, one 
might expect that (( d(t ))) saturates even earlier as re¬ 
pulsive interactions restrict the mobility of the particles 
dynamics even further thereby impeding the emergence 
of any single particle based transport in the system. How¬ 
ever, for a = 2; 5 (( d(t ))) grows linerly, i.e. we observe 
a permanent directed current of the underlying ions that 
emerges spontaneously when the repulsive interaction ex¬ 
ceeds a certain threshold. Fig. [0 inset (b) shows that 
apart from a few ensembles that stay close to their ini¬ 
tial positions, after 1000 oscillations of the lattice (d(t)) 
is very similar for most ensembles which propagate with 
similar velocities through the lattice. For stronger in¬ 
teractions (a = 20; 175) (( d(t ))) saturates quickly as for 
a = 0 , i.e. the directed current dies out after a short 
initial drift and each of the underlying ions is trapped on 
its lattice site. Hence, we observe a directed current that 
is interaction induced but does not survive for strong in¬ 
teractions. 



FIG. 1. Time evolution (( d))(t ) for N = 6 (ochre: N = 20), 
averaged over 10 4 ensembles, for a = 0 (red, 2 nd lowest); 
a = 1 (green, 4th lowest); a = 2 (blue, highest); a = 5 (black, 
2nd highest; ochre 3rd highest); a = 20 (cyan, 3rd lowest); 
a = 175 (purple, lowest). Inset (a): Magnification of the short 
time behaviour of (( d(t))). Inset (b): Distribution of ( d) for 
a = 5.0 and 10000 ensembles at t = 1000(27t/wo). Parame¬ 
ters: a = m = k = 1.0; V = 25; 7 = 8.6; Vo G (—250,250); 
w 0 = 3.53605 (2w 0 = sJWWjm). 



FIG. 2. Asymptotic trajectories Xj(t) (color/grey scale de¬ 
notes the value of Xj(t ) — L/2 modulo L = ir) over particle 
number j = 1,2 ...6 and time t/T. Inset: Low frequency ex¬ 
tract of the squared Fourier spectrum of a typical Vj(t) := 
Xj(t) trajectory. Parameters as in Fig. [T] (N = 6, a = 5). 


It is illustrative to further resolve this effect on a trajec¬ 
tory level (Fig. [2|. While for all non transporting cases 
(a = 0;1;20;175) we asymptotically obtain exclusively 
on site oscillations (not shown), for a = 5.0 the trans¬ 
port is represented by a spatiotemporal pattern. The 
latter consists of phases of on-site oscillations (exemplar- 
ily encircled and marked with an ’O’) that are interrupted 
by propagation processes from one lattice site to its right 
neighbour (’P’) and by propagation processes of two lat¬ 
tice sites (’PP’). The spatiotemporal pattern of the on¬ 
site oscillations and propagation processes is reminescent 
of a situation where all particles oscillate approximately 
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FIG. 3. (a): Asymptotic single particle (a = 0) transport ve¬ 
locity (( Vd)) (blue, upper curve) and kinetic energy per par¬ 
ticle ((-Ekm)} = (m/2)((i 2 )) (green, lower curve) as a func¬ 
tion of 7, respectively averaged over 10 4 single particle en¬ 
sembles and the time interval t/T £ (500,1000). Remaining 
parameters as in Fig. [l] (b): Projection of the asymptotic 
quasiperiodic dynamics of representative, exemplary trajec¬ 
tories in the (x, v, t)-single particle phase space onto the x — v 
plane (x taken modulo L = n) for 7 = 8.0, a = 0 (black, 
highest saturation), 7 = 8 . 6 , a = 0 (blue), 7 = 8 . 6 , a = 5 
(green, lowest saturation), (c): Asymptotic kinetic ener¬ 
gies per particle ((Eki n ))/N for N = 6 and a = 0 (red); 
a = 1 (green); a = 2 (blue); a = 5 (black); a = 20 (cyan); 
a = 175 (purple). The red dot is lying on top of the green 
and the cyan one. (d): Time evolution of the kinetic en¬ 
ergy Ea '■= (m/2)A 2 (t) for a = 0 (red, curves without dip), 
a = 5.0 (blue), A(t = 0) = 2.25 and different A (t = 0) < 0 
(line styles), m = 7 = 1.0. 


in phase while an additional travelling collective exci¬ 
tation represents the propagation process leading to di¬ 
rected transport. 

A Fourier analysis of the transporting pattern (inset 
in Fig. n reveals that it is of quasi-periodic character 
persisting asymptotically. The low frequency part of 
the Fourier spectrum shows subharmonic frequencies of 
nwo/48 with n = 1, 2, 3... and 2ujq = \/2 Vk 2 /m reflecting 
the collective nature of the interaction induced directed 
transport. 

Analysis We now focus on the mechanism underlying 
the interaction induced directed currents. Since our non¬ 
linear system is driven out of equilibrium and breaks both 
parity- and time-reversal symmetry it fulfills the general 
criteria for observing a directed current [s]. Indeed, in the 
single particle regime (a = 0 or N = 1) a directed cur¬ 
rent can be observed, as long as 7 is sufficiently small, i.e. 
for weak dissipation (Fig. [3b; blue curve). This current 
is represented by one or several coexisting transporting 
limit cycle attractors (Fig. [3b; black) in the underlying 
single particle phase space (x,v,t). Physically such at¬ 


tractors exist as long as the energy transfer from the os¬ 
cillating lattice to the particles can periodically balance 
the frictional energy loss, which often happens within a 
single period of the driving. The slowest possible directed 
motion which allows for such a balancing is typically a 
motion of a distance L for a driving period T. In synchro¬ 
nization theory this is a 1:1 (or a k : k, k £ N) phase lock¬ 
ing resonance between the particle motion and the lat¬ 
tice oscillation [39>, 40]. A slower motion of e.g. only L/2 
within T would correspond to a particle motion which 
alternatingly gains and looses energy to the lattice and 
does not allow for a balancing of the energy. Note that 
the existence of such a minimal transport velocity is in 
strong contrast to both the Hamiltonian regime (7 = 0) 
and the regime of strong noise (high temperature/low 
particle mass) where the directed transport emerges as an 
asymmetric (chaotic deterministic or Brownian) diffusion 
which can generally become arbitrarily slow p]. In the 
considered dissipative system, the minimal asymptotic, 
time-averaged transport velocity of one particle in the 
lattice that is typically possible is Vd = 0Jo/(2k) ~ 1.77, 
corresponding to the long plateau in (Fig. [3b; blue) be¬ 
tween 7 ft: 2.1 and 7 ss 7.3. When 7 becomes sufficiently 
large, the transporting attractors (Fig. [3b; black) disap¬ 
pear and we have exclusively non transporting attractors 
(Fig. [3b; blue) in the underlying phase space. Values for 
(( Vd )) that are smaller than u>o/(2k) in Fig. [3b are facil¬ 
itated by non transporting attractors that coexist with 
transporting ones. Some initial conditions then lead to 
on-site oscillations while others lead to transport with 
(vd) = uJo/(2k). Larger transport velocities are possible 
due to higher resonances such as the peak at 7 « 1.1 
representing particles that move 4L within 2T, i.e. a 4:2 
phase locking. 

In contrast to the single particle transport that is based 
on the existence of transporting attractors with a certain 
minimal velocity, the presence of interactions allows for 
the emergence of directed currents as a collective phe¬ 
nomenon. Specifically, while for a = 0 the Fourier spec¬ 
trum consists of peaks at to 0 and its higher harmonics 
2wo; 3 o 7)... (not shown) for a = 5 the Fourier spectrum 
(inset in Fig. [2J shows subharmonic frequencies repre¬ 
senting the possibility of a directed current of the inter¬ 
acting ensemble without requiring that the underlying 
particles propagate by (at least) one lattice site per driv¬ 
ing period. Indeed, in Fig. jT] the transporting ensembles 
propagate with (Vd) = L/(8T) = wo/(16 k) < L/T = 
wo/(2fc). On the level of the 2 N + 1 dimensional N par¬ 
ticle phase space, the corresponding emerging collective 
current (Fig. [l|) is represented by a transporting attrac¬ 
tor (Fig. [3b; green). For visualizing this high dimensional 
attractor we exploit the fact that all N particle trajec¬ 
tories underlying the corresponding transporting pattern 
(Fig. [2]) are identical except for certain phase shifts, i.e. 
the single trajectory shown in Fig. [3b is representative 
for the complete N particle attractor. 
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Let us further illuminate the microscopic mechanisms 
underlying the interaction induced transport. The in¬ 
teraction among the particles can, in the course of the 
dynamics, reduce the loss of kinetic energy due to fric¬ 
tion, which in turn of course facilitates the emergence 
of directed currents. Consider the time evolution of the 
kinetic energy E&(t) := to A 2 / 2 corresponding to the rel¬ 
ative motion A := X 2 — 27 of a free frictious two-particle 
system mx 1,2 + 7x1,2 ± a/[x 1 — X2] 2 = 0. While the fric¬ 
tional energy loss concerning the center of mass dynamics 
E x (t) := mx 2 /2 with x := (27 + X2)/2 is of course inde¬ 
pendent of the value of a, E/±(t) initially decreases faster 
for a > 0 as compared to a = 0 (see Fig. Ol) but then 
grows again and becomes larger than for a = 0 (where 
E/±(t) oc exp [—27^/771]). The zero in E&(t) (see blue 
curves in Fig.jdji) hereby indicates the close encounter ie. 
collision of two particles. Accordingly, a particle travel¬ 
ling from one lattice site to an adjacent one, can reach 
the latter with a higher kinetic energy when encountering 
a collision on its path. Physically this collision induced 
reduction of frictional energy losses is based on the fact 
that upon a collision kinetic energy is converted into in¬ 
teraction energy which is not subject to frictional losses. 
In the case of the complete N particle system in the os¬ 
cillating lattice, this mechanism is still present and leads 
to a substantial increase of the asymptotic kinetic ener¬ 
gies (Fig. [3t) for a = 2; 5 (transport) as compared to the 
energies in the noninteracting a = 0 and non transport¬ 
ing a = 1; 20; 175 cases. In these non transporting cases 
the particles asymptotically exhibit almost equidistant 
on site oscillations. 

We now discuss why the interaction induced directed 
transport is present only in some regime of the inter¬ 
action strength (e.g. for a = 2; a = 5) but vanishes 
for comparatively large values of a (a = 20; a = 175). 
Clearly, in the limit a 00 of a stiff ion chain where all 
N particles exhibit an identical dynamics, the N parti¬ 
cle transport breaks down at the same friction strength 
7 as in the single particle case (N = 1 or a = 0). Us¬ 
ing linear stability analysis, we show that the N particle 
dynamics reduces to that of a stiff ion chain already for 
finite values of a, no matter how strong the driving is. 
Consider a two particle system (Eq. ([]]) for N = 2; r = 1) 
with periodic boundary conditions and only next neigh¬ 
bour interactions. Transforming 27 —> 27 := 2fc[27 — /(f)] 
and then to center of mass x := (27 + X 2)/2 and relative 
A := X 2 — X 1 coordinates yields after inserting A = 27r+<5 
(2-7T is the equilibrium distance) and truncating in linear 
order 5: 


0 = mx + 72; — 2Vofc 2 sin(x) + F(t) 

0 = md -t-7<5 + 2Vok 2 6 cos(2x) + 8k 3 a6/n 3 


F(t) = 2k mf(t ) + 7 f(t) 


( 2 ) 

(3) 


As cos(2x) £ (—1,1) it follows from the linear Eq. J3]) 
that for a > V^tt 3 /(Ak) the asymptotic state for S(t) is 



FIG. 4. Time evolution (( d))(t ) for N = 6 , averaged over 10 4 
ensembles for a = 0.25; 2; 5 (three upper curves in green, blue 
and black, lying almost on top of each other); a = 20 (cyan, 
second lowest curve); a = 175 (purple, lowest curve) and 
a = 0 (covered by the purple curve). Inset (a): Like Fig. [2] 
Inset (b): Like Fig. [3J;.) Upper dots: a = 0.25 (covered by 
the blue dot), a = 2 (blue) and a = 5 (black). Lower dots: 
a = 20 (cyan) and a = 175 (purple, covering the a = 0 dot 
and partly also the cyan one) 

Parameters: ujq = 4.0; a = m = 1.0; 7 = 7.0; V = 25; vo G 
(-250,250). 


uniquely determined by a fixed point attractor at S = 0 
with a global basin of attraction. Then, the state where 
the two considered ions exhibit an identical phase locked 
dynamics is stable (for arbitrary values of a, to 0 and 7) 
and any initial state sufficiently close to it asymptoti¬ 
cally approaches it. This result can be generalized to 
the N particle case, yielding a > KVqtt 3 /(Qk) with 
Q = 4; 3; 2; (5 - v / 5)/2; 1 for N = 2,3; ...6 and Q > 0 
for N < 00. In fact, for a > KVqtt 3 /(Qk) we numeri¬ 
cally observ that the N particle ensemble discussed in the 
context of Fig. |T]generically comes close to a translation- 
ally invariant state (after redistributing the nonuniform 
initial kinetic energies) and we observe a breakdown of 
the N particle current at the same value of 7 as in the 
single particle case. The uniform (or phase locked) N 
particle dynamics is then described by the equation of 
motion of its center of mass coordinate x := (1/iV) 27 

(Eq.[2]), representing a biharmonic version of an ac-driven 
damped physical pendulum which is known to exhibit a 
period doubling route to chaos [+lJ. Note that the cri¬ 
terion a > KVqtt 3 /(Qk) represents only an upper limit 
for the stability of the spatially uniform dynamics; in 
practice stability can of course occur already for smaller 
values of a. 

Irregular transport We finally demonstrate that the 
interaction induced directed currents are by no means 
restricted to the regular regime as represented by the 
propagating spatiotemporal pattern in Fig. [21 Fig. [4] 
shows (( d(t ))) in a regime where the dynamics of the 
underlying ions is irregular (Fig. [4l inset a). Here, a 
current emerges already for comparatively weak interac¬ 
tions (a = 0.25) and is accompanied, as before, by a 
significant increase of the mean kinetic energy per parti¬ 
cle compared to the noninteracting (a = 0) and the other 
nontransporting cases (a = 20; a = 175) (Fig. [2 inset 
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b). The corresponding Fourier spectrum of the velocities 
of the ions (not shown) exhibits a continuous distribu¬ 
tion. This provides the particle system with an enhanced 
flexibility to dynamically transfer energy between the in¬ 
dividual particles. The emerging current is larger than 
that of the transporting regular excitation pattern dis¬ 
cussed {((d)) ~ 0.2wo/fc) but still smaller than any sin¬ 
gle particle transport corresponding to a 1:1 resonance 
((d)) = u] 0 /(2k). 


Conclusions Opposite to the common assumption 
that directed currents represent a pure single particle ef¬ 
fect that can survive interactions, we demonstrate here, 
that collective effects in repulsively power-law-interacting 
systems can induce directed particle currents even in pa¬ 
rameter regimes where corresponding noninteracting sys¬ 
tems do not exhibit transport. These interaction induced 
currents could be detected e.g. for cold thermal clouds 
of dipolar atoms and molecules 3C)|, 31 and interacting 
colloids 32. •A- in optical lattices driven via standard 
techniques 31 M or alternatively with trapped ions ad¬ 
ditionally exposed to driven optical lattices [33]. Our 
results might open a doorway towards the largly unex¬ 
plored territory of interaction based transport phenom¬ 
ena in time dependent lattices. 
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